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ABSTRACT 


A  theoretical  analysis  of  the  axisymmetric  elastic  deformations  and 
stresses  in  a  web*stiffened  sandwich  cylindrical  shell  structure  under  ex* 
ternal  hydrostatic  pressure  is  presented.  The  solution  is  based  on  the  use 
of  edge  coefficients  for  plate  and  shell  elements  of  finite  length,  and  in¬ 
cludes  the  computation  of  the  edge  forces  and  moments  arising  at  the 
common  junctures  of  these  elements. 

Equations  are  given  for  computing  numerically  the  longitudinal  and 
circumferential  stresses  in  the  two  coaxial  cylindrical  shells  and  the  radial 
and  tangential  stresses  in  the  web  stiffeners  between  the  two  shells. 

No  consideration  was  given  to  the  discontinuity  effects  arising  from 
rigid  or  elastic  restraints  afforded  by  contiguous  bulkhead  or  adjacent  shell 
structures.  Thus,  the  analysis  presented  herein  is  applicable  only  to  a 
typical  bay  of  a  web-stiffened  sandwich  cylinder  of  long  length. 

A  numerical  example  is  presented  to  illustrate  the  use  of  the  equations 
developed  in  this  report. 


INTRODUCTION 

The  David  Taylor  Model  Basin,  under  initial  sponsorship  by  the  Office  of  Naval 
Research  and  later  continuance  by  the  Bureau  of  Ships,  has  been  investigating  the  feasi¬ 
bility  of  sandwich-type  construction  for  pressure  hull  application.  Results  of  exploratory 
experimental  studies  carried  out  under  this  program^  have  shown  that  in  certain  ranges  of 
geometry  strength-weight  advantages  on  the  order  of  20  to  25  percent  higher  can  be  realized 
with  sandwich  designs  over  the  conventional  ring-stiffened  cylindrical  configuration.  These 
results  were  obtained  from  model  tests  of  sandwich-type  cylinders  having  "hard”  cores;  i.e. 
the  cores  were  capable  of  developing  high  compressive  strengths  in  addition  to  transmitting 
the  pressure  loading  by  shear  from  the  outer  to  the  inner  shell. 

At  the  time  these  sandwich  cylinders  were  conceived,  no  formulas  were  available  on 
which  to  base  an  optimum  design;  merely  intuition  and  engineering  judgment  were  resorted 
to  for  proportioning  the  structural  elements.  Concurrently  with  the  experimental  program, 
analytical  studies  were  initiated  to  develop  rational  formulas  based  on  thin-shell  theory  for 
predicting  the  elastic  deformations  and  stresses  in  the  structural  elements  of  such  sandwich 
type  cylinders. 

In  this  report,  equations  are  developed  for  carrying  out  a  complete  stress  analysis  of 
a  typical  portion  of  a  web-stiffened  sandwich  cylinder  under  external  hydrostatic  pressure. 

^Refarencas  are  listed  on  page  37. 


The  method  is  based  on  the  use  of  edge  coefficients  of  plate  and  shell  elements  of  finite 
length,  and  satisfaction  of  force  and  moment  equilibrium  and  compatibility  of  deformations 
at  the  common  junctures  of  the  elements  comprising  the  structure.  Expressions  for  edge 
coefficients  of  cylindrical  shells  of  short  length  are  developed  in  Appendix  A. 

GENERAL  CONSIDERATIONS 

Methods  of  analysis  based  on  the  use  of  edge  coefficients  have  found  wide  application 
in  studying  stresses  and  deformations  in  complex  structures  composed  of  ring,  plate,  and  shell 
elements. The  underlying  concept  in  this  type  of  analysis  is  that  a  complex  physical 
structure  can  be  broken  down  into  identifiable  components  for  which  mathematical  solutions 
exist  or  can  be  found  readily.  The  deformations  occurring  in  each  structural  element  are  deter¬ 
mined  in  terms  of  unknown  forces  and  moments  assumed  to  exist  at  the  junctures  common  to 
these  elements.  Conditions  of  equilibrium  and  compatibility  are  then  satisfied  at  each  of  the 
junctures,  thus  permitting  determination  of  the  redundant  forces  and  moments.  With  this 
information,  a  complete  stress  analysis  for  each  structural  component  can  then  be  performed. 

The  present  problem  of  the  stresses  in  a  web-stiffened  sandwich  cylinder  subjected  to 
hydrostatic  pressure,  shown  in  Figure  1,  can  be  solved  rather  conveniently  by  the  use  of  edge 
coefficients.  The  identifiable  structural  elements  in  this  case  are  two  coaxial  cylindrical 


Figure  1  -  Web-Stiffened  Sandwich  Cylinder 
Subjected  to  External  Hydrostatic  Pressure 


shells,  one  subjected  to  radial  pressure  and  an  end  load  and  the  other  to  an  end  load  only, 

a 

and  annular  discs  subjected  to  radial  loads  on  the  two  circular  boundaries.  The  webs  or 
annular  discs  act  as  the  connecting  and  stiffening  members  to  the  two  shells.  A  free-body 
diagram  showing  the  breakdown  of  the  physical  structure  to  its  component  parts,  together 
with  appropriate,  but  as  yet  unknown,  edge  forces  and  moments,  is  presented  in  Figure  2. 


2 


Figure  2  -  Free-Body  Diagram  Showing 
Forces  and  Moments  Acting  on  Shell 
and  Web  Elements  of  Web*Stiffened 
Sandwich  Cylinder 


Following  the  method  of  References  2,  3,  and  4,  the  deformations  occurring  at  the 
edges  of  a  shell  element  of  general  meridional  shape  can,  by  simple  superposition,  be  written 
in  terms  of  the  unknown  edge  forces  and  edge  moments  and  known  applied  loading  as  follows: 


^  flP,  4  /;>,  4/;"P.  4  d,.M.  4  g,.Q.  4 
d.  =  a.M.  4  6.(3,.  4  cfp.  4  Cf'p  4  cpp.  4  a.jM.  4  6, ..(3^.  4  c,rP. 


[1] 

[2] 


where  the  coefficients  a,.,  6,.,  .  .  .  are  the  amount  of  transverse  deflection  or  meridional 
rotation,  as  the  case  may  be,  per  unit  bending  moment,  shearing  force,  axial  force,  or  surface 
pressure  loading,  as  shown  in  Figure  3.  The  coefficients  with  the  double  subscripts,  i.e.. 


Figure  3  -  Shell  Element  of  Arbitrary 
Meridional  Shape  Subjected  to  Edge 
Moments,  Shears,  Forces,  and 
Surface  Loading 


a.j,  b-p  .  .  ,  flj,  are  the  interaction  coefficients  which  reflect  the  deformations  at  edge  “z" 
due  to  forces  and  moments  at  edge  By  replacing  i-*j  and  }-»i  in  Equations  [1]  and  [2], 

expressions  for  the  deformations  Wj  and  6-  can  be  written  immediately. 

Note  that  the  effect  of  the  end  load  P  on  the  deformations  u),.  and  6-  has  been  separated 
into  three  distinct  components.  The  components  denoted  by  the  single-primed  coefficients 
/,^and  c,r  are  those  due  to  bending  effects.  The  same  is  true  of  the  components  associated 
with  the  coefficients  /,w  and  c,y,  but  these  also  reflect  interaction  influences.  The  components 


3 


denoted  by  the  triple-primed  coefficients  /f"and  c-"&re  essentially  Poisson  effects  on  the 
membrane  deformations. 

For  the  specific  problem  of  cylindrical  shell  elements  symmetrically  loaded,  as  shown 
in  Figure  4  and  considered  in  this  report,  Equations  [1]  and  [2]  become: 


M. 


+  P 


dMx 

dx 


-l-Mx  ,  f  M.  s  M 

2  Qx  «  H| 
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|/2 


Figure  4  —  Sign  Convention  for  Cylindrical-Shell  Element  (Symmetric  Case) 


et  =  ^  clP^  crff  ^  ^  c;.p, 


[3] 

[4] 


and  similar  expressions  for  and  O.,  respectively.  However,  for  the  case  of  a  cylinder  some 
of  the  terms  appearing  in  Equations  [3]  and  [4]  become  zero;  this  will  be  shown  later.  In 
addition,  for  the  pressure  loading  shown  in  Figure  1,  where  the  inner  cylindrical  shell  is  not 
subjected  to  the  radial  pressure  loading,  those  terms  in  Equations  [3]  and  [4]  that  are  multi¬ 
plied  by  p  will  drop  out  when  the  deformations  of  the  inner  shell,  i.e.  i  =  2,  are  considered. 

Following  the  same  technique  employed  for  the  shell  elements,  the  deformations 
occurring  at  the  edges  of  the  circular  annuli  or  web  elements,  as  shown  in  Figure  5, 


I 


Figure  5  —  Sign  Convention  for 
Web  Element 
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can  be  written  in  terms  of  the  unknown  edge  forces  and  moments  and  known  applied  loading 
as  follows: 


[5] 

of  +dfp  +  [6] 

where  it  has  been  assumed  that  the  edge  moments  Mf ,  axial  thrusts  and  surface  pressure 
p  do  not  give  rise  to  any  radial  displacement  w  in  the  plane  of  the  circular  annulus.  Equations 

[5]  and  [6]  give  the  deformations  at  boundary  i  of  the  circular  annulus;  the  deformations  at  the 
other  boundary,  say  can  be  obtained  by  replacing  i  by  j  and  ;  by  i  in  Equations  [5]  and  [6], 
The  expression  for  the  edge  rotations  df  and  of  the  annulus  are  rather  general  to  include 
the  case  in  which  the  sandwich  void  between  the  two  cylinders  may  become  pressurized.  This 
problem  will  not  be  considered  in  this  report.  In  a  later  section,  it  will  be  shown  that  due  to 
symmetry  the  edge  rotations  of  the  annulus  are  zero.  Furthermore,  for  the  particular  case  of 
pressure  loading  shown  in  Figure  1  and  considered  in  detail  later  including  a  numerical  ex* 
ample,  not  only  are  the  edge  rotations  df  and  0^  equal  to  zero  but  every  term  in  expressions 

[6]  is  zero.  In  such  a  case  it  is  tacitly  assumed  that  the  web  stiffeners  act  only  to  resist 
hoop  compression  and  do  not  act  in  the  sense  of  a  circular  plate  to  resist  bending  due  to  edge 
moments  and  edge  shears. 

In  Reference  4,  for  instance,  equations  were  developed  for  computing  discontinuity 
stresses  at  cone-cylinder  junctures,  either  with  or  without  transverse  reinforcement.  For 
that  problem  it  was  tacitly  assumed  that  the  shell  elements  were  each  of  semi-infinite  length 
so  that  the  deformations  at  their  common  juncture  were  not  influenced  by  boundary  effects  at 
the  others  ends.  This  permitted  the  use  of  rather  simple  expressions  for  the  edge  deformations. 

For  the  present  problem  of  the  web-stiffened  sandwich  cylinder,  the  elements  com¬ 
prising  the  structure  are  of  such  proportions  that  interaction  of  internal  edge  effects  is  very 
predominant.  This  necessitated  the  development  of  edge  coefficients  for  cylindrical  shell 
elements  of  finite  length.  However,  it  turns  out  that  the  forms  of  the  new  coefficients  are 
exactly  the  same  as  those  of  Reference  4  except  for  multiplying  factors  which  are  functions 
of  the  shell  geometry  and,  primarily,  the  length.  These  edge  coefficients  for  a  cylindrical 
shell  are  written  in  the  following  convenient  form: 


£*0  =+  -j- 


£*6.= - — 


E^c;=  Ecl'=Ecl"=  Ef;=0 
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1 


=  + 


2£?;/3f 


Evr=+ 


Rf 

h, 


R, 

E^fr--^  — 

fli 


E^g.  = - L_  At3](e.O 

2d;b^ 


[7a] 


£®a..  - - a[.®\/3.7) 

*/  0/(8,.  ‘  ‘ 

EH..  =  +  — - — 

^  ^n'o2 


2D;^f 


E^c,.^Eff.  =  0 


2Dl^} 


where 


A’  12(1-1,2) 


!L_  .  fl  ^3(1 -^,2) 


[7b] 


The  “lambda”  functions  A^^^,  a[^^,  a[.^\  .  .  .  A^®^  appearing  in  the  edge  coefficients, 
Equations  [7a],  are  derived  in  Appendix  A  and  are  defined  here  as  follows: 


aS'\/3.0  = 


sinh^  ^.l  +  sin^  j8,i 


AS'^(i3,0  = 


sinh^ -  sin^  /3,f 

cosh  sinh  +  cos  sin  j8,l 
sinh^/3,f  -  sin^/Sj.Z 
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cosh sinh  -  cos  sin 

a!>'(s,,o - - - ^ - 

sinh^ p.l  - 

[8] 

^  ,  2sinh Bil  sin  i3,Z 

- ^ 

sinh^  6fl  -  sin^  B^l 

cosh  sin /3,Z  -  sinh /3,7  cos 

AS>l((i,i) - - - j - 

sinh^  -  sin^ 

cosh ^.l  sin  ^.l  +  sinh j3  I  cos  /S  I 

aI»1o,!) - — - J - 

sinh^  /Sj.Z  -  sin^  B^l 


For  convenience  and  ease  of  calculation,  numerical  values  of  the  “lambda”  functions  in  [8] 
were  determined  with  the  aid  of  a  Burroughs  E-101  computer  for  a  range  of  13 -I  from  0.40  to 
2.50  in  increments  of  0.02.  The  results  were  tabulated  and  are  given  in  this  report  as  Table  1. 

For  the  special  case  of  a  cylindrical  shell  of  semi-infinite  length,  i.e.,  ?-»«>,  the 
interaction  functions  given  by  Equations  [8]  simplify  to 


a[i]  =  a[2]  =  aW  =  1 
=  a[®^  =  a[.®^  =  0 


[9] 


and  the  edge  coefficients  given  by  Equations  [7a]  reduce  exactly  to  those  given  in  Reference  4. 

From  symmetry  considerations  it  is  seen  that  the  edges  of  the  web  stiffener,  which  for 
purposes  of  analysis  is  viewed  as  a  circular  annulus,  do  not  undergo  any  rotation.  This  stems 
from  the  fact  that  a  horizontal  tangent  or  zero-slope  condition  is  assumed  to  exist  at  the 
junctures  of  the  webs  with  the  two  cylindrical  shells.  This  assumption  implies  that  the  edge 
moments  on  each  shell  at  the  shell-web  junctures  balance  each  other,  so  that  there  are  no  net 
moments  to  be  resisted  by  the  web.  Further,  it  is  assumed  that  the  web  elements  do  not  take 
any  axial  force  due  to  the  axial  pressure,  but  that  this  is  all  resisted  by  the  cylindrical  shells. 
Thus,  the  analysis  of  the  web  stiffener  is  reduced  to  that  of  a  circular  annulus  subjected  to 
axisymmetric  in-plane  radial  forces  on  both  its  inner  and  outer  boundaries;^  see  Figure  2. 

On  the  basis  of  these  assumptions,  it  is  necessary  to  derive  edge  coefficients  for  an  annulus 
undergoing  radial  deflections  only.  Such  coefficients  are  developed  in  Appendix  B  and  are 
given  here  as  follows: 
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UO 

CM  CM  CM  CM  CM  CM 


where 


(1  +  v)  /?2 

g  a  - — - — - 

A3(/?2-/?2) 


[10a] 


[10b] 


COMPUTATION  OF  STRESSES 

The  formulas  given  herein  for  determining  the  longitudinal  and  circumferential  stresses 
in  the  shell  elements  of  the  web-stiffened  sandwich  cylinder  are  developed  in  Appendix  C. 
Formulas  for  the  radial  and  tangential  stresses  in  the  web  elements  are  developed  in  Appendix 
B.  The  derivation  follows  very  closely  the  general  analysis  of  Reference  7  for  ring-stiffened 
cylinders  under  hydrostatic  pressure,  the  only  differences  arising  from  the  elastic  restraints 
at  the  shell  edges  and  the  distribution  of  the  axial  pressure  loading. 

The  nomenclature  and  sign  convention  used  in  Reference  7  and  in  Appendix  C  of  this 
report  are  shown  in  Figures  2  and  4.  A  longitudinal  bending  moment  is  considered  posi¬ 
tive  if  it  tends  to  put  the  outer  surface  of  the  shell  in  tension,  and  a  transverse  shearing  force 
is  considered  positive  when  it  acts  in  a  direction  away  from  the  axis  of  symmetry  but  in 
the  positive  avdirection.  A  hydrostatic  pressure  p  is  considered  positive  when  it  is  external, 
and  negative  when  internal.  With  reference  to  Equations  [1]  and  [2],  the  subscript  i  is  used 
to  distinguish  the  two  cylinder  elements. 

The  quantities  H-  and  shown  in  Figure  2  are  the  edge  shearing  forces  and  bending 
moments  arising  at  the  junctures  of  the  shell  elements  with  a  web  stiffener.  They  may  be 
determined  in  terms  of  the  geometry  and  elasticity  of  the  structure  and  the  pressure  loading 
by  enforcing  conditions  of  force  and  moment  equilibrium  and  compatibility  of  deformations 
at  the  junctures.  This  determination  is  developed  in  the  next  section. 
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Once  the  edge  forces  and  moments  are  known,  the  following  formulas  may  be  used  for 
determining  the  critical  longitudinal  and  circumferential  shell  stresses  which  occur  at  a  point 
midbay  between  two  adjacent  webs  and  also  at  a  web  location,  respectively: 


AT  MIDBAY: 


^Xm 


(-«*) 


£:*A../32  A[®\i8,7/2) 

-  (tflP-u)?)  — — - 

— - + - - 


a\^\^,1/2) 


[11] 

[12] 


AT  A  WEB: 


^Xf 


P.  a[^\^.1/2) 

—  ±  -  K-O  - 


A; 


(1-v") 


= 


A, 


«.• 


uE»h.^f 

(1-^2) 


a\^\^,1/2) 

At3](^.f/2) 


(wP-w?) 


A\^\^il/2) 


[13] 

[14] 


where  in  the  above  equations  j  «=  1,  2,  and  the  upper  sign  is  for  the  outer  fiber  and  the  lower 
sign  for  the  inner  fiber  of  each  shell  plating.  Equations  [11]  through  [14]  are  developed  in 
Appendix  C. 

Once  the  critical  stresses  are  determined  from  Equations  [11]  through  [14],  the 
question  as  to  how  they  combine  to  precipitate  axisymmetric  collapse  of  the  cylindrical  shell 
elements  can  be  answered  by  recourse  to  the  failure  criteria  discussed  in  References  7  and  8. 
This  will  not  be  discussed  here. 

The  quantity  P.  in  Equations  [11]  through  [14]  is  the  axial  load  taken  by  each  of  the 
two  cylindrical  shells.  On  the  assumption  that  the  two  shells  contract  the  same  amount 
longitudinally,  it  is  shown  in  Appendix  D  that 
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uE^h.  ^  vE*h^ 

pf^l 

D  1 

+  p  (^1  ^2)^1“  p  ^^2  ^2*^2  ^1)^2 

it  j  ^2 

“ 

R.h. 

1 

[15] 


p«i 


^2- 


r  R,  1  R.  ^ 

-.(1-12,)  Y*'iT  ('-X)  -— 

2  2  I  ^ 


«? 


1  -  v^fl,  + 


^  ^2^2 


(l-v2fl,) 


The  quantities  <?j  and  are  given  in  Appendix  D  by  Equations  [D.13]. 

The  quantities  w?  in  Equations  [11]  through  [14]  represent  the  particular  integrals  to  the 
differential  equations  governing  the  axisymmetric  deformations  of  a  cylindrical  shell;  they  are 
easily  determined  from  membrane  theory.  For  the  case  shown  in  Figures  1  and  2,  where  the 
outer  cylindrical  shell  (but  not  the  inner  one)  is  loaded  by  lateral  pressure,  we  find  that  (see 
Reference  7,  for  example) 


VRI 

(1-  ^ 

"A 

Wf  ^  +  ! 

E*Aj 

\ 

V  / 

2,  1 

- — 

E’hj 

f  V  ^2 

\ 

\^2~P 

I 

where  the  axial  forces  and  Pj  are  given  by  Equations  [15], 

The  shell  edge  deflections  w*  appearing  in  the  stress  formulas,  Equations  [11]  through 
[14],  are  determined  from  Equation  [3]  once  the  edge  shears  and  edge  moments  are 
known;  i.e., 

- 

^  ^  ^  ^  ^  T  ^ 

[17] 

H  H 

tOj  =  <^2^2  +  ^2  —  +  f2^2  +  ^2^2  ^2  ^2^2 

where  the  interaction  coefficients  have  been  designated  by  a  “bar”  instead  of  the  double  sub¬ 
script  so  as  not  to  confuse  the  use  of  the  subscripts  “1”  and  “2”  to  designate  the  two  shells 
and  their  respective  junctures  with  the  web  stiffeners.  This  notation  will  be  used  in  all  the 
equations  that  follow. 
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Expressions  for  the  radial  and  tangential  stresses  in  the  web  elements  are  developed 
in  Appendix  B.  It  is  shown  there  that  the  maximum  radial  stress  occurs  at  the  intersection 


maximum  tangential 


/  -  ^\ 

with  the  outer  cylindrical  shell  I  i.e.,  at  r  =  /?,  =  ft,  +  —  I  ,  and  the 

\  2  /  /  _  rtj  V 

stress  occurs  at  the  intersection  with  the  inner  cylindrical  shell  ^  i.e.,  at  r  »  ftj  -  ftj  — —  I  ■ 


These  maximum  stresses  are  given  by  the  following  expressions: 


(1 


where  the  constants  A  and  B  are  given  by 


E^  [  fi  "1 

-  A(l  +  v)-  — (1-v) 

-u^)  L  J 

E^  V  fi  "1 

— _  (i-v) 

L  Rl  J 


\  ft?  -  ft ?  / 


=  -«l«2  [— - = - 

\  ft?  -  ft?  / 


fi 


[18] 


[19] 


[20] 


and  the  annulus  edge  deflections  and  by  Equation  [5]  as 


w 


[21] 


In  Equations  [21]  the  edge  coefficients  designated  by  a  “bar”  are  the  interaction  or  double¬ 
subscript  coefficients;  i.e.,  s  g^2  ®nd  g^  =  g^i'  The  edge  coefficients  appearing  in 
Equations  [21]  are  given  by  Equations  [10]. 


DETERMINATION  OF  EDGE  SHEARS  ft,.  AND  EDGE  MOMENTS  ft,. 

For  the  case  of  symmetry  on  each  side  of  a  web  stiffener,  the  conditions  of  force  and 
moment  equilibrium  at  each  of  the  two  junctures  of  the  web  with  the  shells  are  rather  obvious; 
these  are  shown  in  the  free-body  diagram  of  Figure  2.  There  remains  to  determine  the  unknown 
edge  shears  ftj  and  and  unknown  edge  moments  ftj  and  ftj  by  enforcing  conditions  of  com¬ 
patibility  of  the  deformations  at  the  junctures  labeled  “1”  and  “2.” 
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Continuity  and  symmetry  conditions  at  joint  “1”  require  that 


1C*  = 

61  =  6^  ^0 

whereas  these  conditions  applied  to  joint  “2”  require  that 

icj  = 

«*  =  6^  =0 


[22] 

[23] 


[24] 

[25] 


Substituting  Equations  [3],  [4],  [5],  and  [6]  into  the  four  conditions  [22]  through  [25], 
considering  the  zero  edge  coefficients  by  virtue  of  the  loading  shown  in  Figures  1  and  2  and 
Equations  [7],  and  assuming  that  £*  i:  S'*,  we  must  solve  the  following  four  algebraic  equations 
simultaneously  to  determine  and 


Af  -  gf^  +  -  -  /,'>  -  f{"P,  [26] 

[aj  +  Oj]  +  Wj  (6j  +  6j)j  =  0  [27] 

[d^^d^]  +  (g^^g^)  -  g^  +  [-g^]  =  -  [28] 

^2  [oj+oj] + //j  ® 


where  and  P^  are  given  by  Equations  [15]  to  be  functions  of  the  unknown  shearing  forces 
//j  and  H2’  Equations  [26]  through  [29]  can  be  rewritten  as  two  equations  in  only  two 
unknowns  as  follows: 


«i \j  -  i'.'’  -  J  -  «2  l-f'i'I  -  -  /.->  -  /.'"'1 1301 


"■[i 


H^l-g^]  +  H2  T  (92*9^-92  - 


(i 


5  +  62)  (‘^2+‘^2)"| 

2(02  +  02)  J 


f2"P2  [311 
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By  the  substitution  of  the  expressions  for  Pj  and  given  by  Equations  [15]  into  Equations 
[30]  and  [31],  the  two  simultaneous  equations  to  be  solved  for  and  Pj  become: 

«!  _  \ 


H, 


(6^+ 6j)(dj +  dj)  vE^h^ 

'  —  ■  +  ' 
2  (Oj  +  flj)  R^ 


f: 


‘  «2^2 


(1-v 


^^2)J 


P, 


■ffi  + 


fx 


R.h. 


-n'r-f'x 


r  1  *  A  1 


^2^2 


(l-v2fl-) 


[32] 


Hr 


H, 


-9^  - 


vP®A, 


P  2 
"2 


P, 


P. 


^292 


PiAj 


^  PjAj 


(i-v2a,)j 


/  +  -:r- 


2  (i72-^2)-^2^-  2(0, .a-,) 


Rn 


f'" 

'2 


•W 


P. 


^292 


=  -/2 


r  1  1 

^  J  ro  oi 


R,h, 


^  P2P2 


After  the  edge  shears  and  P,  are  determined  from  Equations  [32]  and  [33],  the  edge  moments 
and  P,  may  be  found  from  the  following  expressions  as  a  consequence  of  Equations  [27]  and 
[29],  respectively: 

(^  +  ^) 

Pi  =  _  Pi  - -  [34] 

2  (Oj  +  Oj) 

(62+62) 

M2  =  ~  P2  — .  [36] 

2(02  +  0,) 
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NUMERICAL  EXAMPLE 


As  part  of  its  research  and  evaluation  program  to  study  the  application  of  glass-fiber 
reinforced  plastics  for  pressure  vessel  construction,  the  Model  Basin  in  collaboration  with 
Narmco  Industries  Inc.,  San  Diego,  California,  is  presently  designing  a  series  of  web-stiffened 
sandwich  cylinders  made  of  these  materials.  The  structural  models  are  to  be  fabricated  by 
Narmco  and  then  forwarded  to  the  Model  Basin  for  testing. 

One  of  the  designs.  Model  N-1,  presently  being  conceived  will  be  used  as  a  sample 
calculation  to  illustrate  the  use  of  the  equations  developed  in  this  report.  The  detailed 
dimensions  are  shown  in  Figure  6,  and  are  summarized  here: 


Figure  6  -  Schematic  Diagram  Showing 
Dimensions  of  Model  N-1 


Aj  »  0.142  in.;  Aj  0>124  in.;  Aj- 0.107  in. 
=  3.493  in.;  =*  3.062  in.;  I  =  0.648  in. 

=  3.564  in.;  »  3.000  in. 

E*  -  6.0  X  10®  psi  - 
V  -  0.15 


Using  Equations  [7b]  for  each  of  the  outer  and  inner  cylindrical  shells,  respectively, 
we  compute  the  values  of  and  to  be; 


=  1.204;  £>j'-  24.407  X  10"®  in.3 
/Sjf  -  1.376;  D'~  16.257  x  10"®  in.® 


The  lambda  functions  are  either  computed  by  using  these  values  of  and  Equations  [8]  or 
are  found  by  interpolation  from  Table  1  for  each  of  the  two  shells.  They  are  summarized  here: 


aW  a[®^ 

Outer  Shell  2.219  2.162  1.694  1.981  0.8061  1.567 

Inner  Shell  1.778  1.655  1.502  1.471  0.6906  0.9806 
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Next,  the  shell-edge  coefficients  E^a^,  E‘bf,  E*c^,  .  .  .  etc.,  are  computed  by  using 
Equations  [7a].  The  numerical  values  thus  found  are  summarized  here: 


Shell  Edge 

Outer  Shell 

Inner  Shell 

Coefficient 

(i  -  1) 

(<-2) 

£*0,.,  in.“^ 

+  0.4767  X  10“ 

+  0.4795  X  10“ 

£*6j,  in."  ^ 

-  0.1317  X  10“ 

-0.1213  X  10“ 

£*dp  in."^ 

0 

0 

+  0.1317  X  10“ 

+  0.1213  X  10“ 

E*f{ 

0 

0 

E>fr,  in. 

+  85.9229 

nonexistent 

E*f.'" 

-  3.6898 

-3.7040 

-  0.05410  X  10“ 

-0.04822  X  10“ 

E^Ep  in."^ 

-  0.3455  X  10“ 

-0.2840  X  10“ 

E*b.,  in."^ 

+  0.1176  xl0“ 

+  0.1003  X  10“ 

£*^.,  in."‘ 

-  0.1176  X  10“ 

-0.1003  X  10“ 

E»f{ 

0 

0 

E‘gi 

+  0.02575  X  10“ 

+  0.02217  X  10“ 

The  web  stiffener  or  circular  annulus  edge  coefficients  E'^g^ ,  .  .  .  etc.,  are  computed  by 
using  Equations  [10],  and  the  numerical  results  found  are: 

E^g^  ~  190.261 
E'^g^  m  E'^g^^  -  161.949 
E‘*g^  -  168.564 
E‘*g^  m  E^g^^  -  192.396 

The  components  of  the  end  pressure  loading  taken  by  each  of  the  outer  and  inner 
cylindrical  shells,  respectively,  are  computed  to  be,  using  Equations  [15]: 

Pj  -  +  0.9843p  Ib/in. 

Pj  -  +  0.8696 p  Ib/in. 

With  all  this,  the  edge  shear  forces  and  Pj  computed  by  solving  Equations  [32] 
and  [33]  simultaneously.  The  values  thus  found  are  then  substituted  into  Equations  [34]  and 
[35]  to  determine  the  edge  bending  moments  and  Pj'  numerical  values  thus  found  are: 

Pj  -  +  0.3893 p  Ib/in. 

Pj  -- 0.2717  p  Ib/in. 
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/Wj  -  +  0.02095  p  in.-lb/in. 

~  0.01459  p  in.-lb/in. 

When  the  edge  shear  forces  and  edge  bending  moments  are  known,  the  edge  deflections 
of  the  two  cylindrical  shells  and  those  of  the  web  stiffener,  at  their  common  juncture  points, 
are  found  from  Equations  [17]  and  [21]  to  be: 

E*w*  =  +  30.072  p  Ib/in. 

E^Wy  “  +  30.072  p  Ib/in. 

E*w*  =+  29. 105 p  Ib/in. 

E^w^  ■  +  29.105  p  Ib/in. 

Comparison  of  E*Wy  with  E^Wy  ^  and  with  E^w^  affords  a  check  on  the  numerical 

calculations,  since  the  boundary  conditions  [22]  and  [24]  enforced  at  the  two  junctures 
require  them  to  be  equal  in  their  respective  cases  since  it  was  assumed  that  E^  =  E^. 

The  maximum  radial  and  tangential  stresses  in  the  web  stiffeners  can  now  be  com¬ 
puted  by  using  Equations  [18],  [19],  and  [20].  The  values  found  are: 

a,  --  3.638  p  Ib/in.  2 

Oy  --  10.083  p  Ib/in.  2 
max 

Before  the  shell  stresses  can  be  computed,  it  is  necessary  to  determine  the  membrane 
deflections  of  the  two  shells.  This  is  done  with  the  aid  of  Equations  [16].  The  values  found 
are: 

■  +  82.291  p  Ib/in. 

3.2205 p  Ib/in. 

Finally,  the  critical  longitudinal  and  circumferential  shell  stresses  at  points  midbay 
between  two  adjacent  web  stiffeners  and  at  a  web  stiffener  are  determined  by  using  Equations 
[11],  [12],  [13],  and  [14].  The  numerical  values  are  summarized  as  follows: 


psi 

psi 

Outer  Shell 

-10.039  p 

- 10.437 p 

-13.166p 

- 10.584 p 

Inner  Shell 

-  9.843p 

- 12.702  p 

-12. 704  p 

- 13.514  p 

For  the  numerical  example  considered,  the  calculations  already  carried  out  have  been 
based  on  the  assumption  that  all  structural  elements  have  the  same  elastic  modulus  E.  How¬ 
ever,  in  the  fabrication  of  a  shell  structure  such  as  this,  it  is  conceivable  that  the  elements 
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could  have  different  material  properties.  In  the  case  of  Model  N-1,  which  is  to  be  made  of  a 
glass-fiber  reinforced  plastic,  it  is  expected  that  the  web  stiffeners,  although  made  of  the 
same  basic  material  as  the  cylindrical  shells,  will  have  a  higher  elastic  modulus  by  virtue  of 
the  fiber  distribution.  Assuming  that  the  modulus  of  the  web  material  is  50  percent  higher 
than  that  of  the  shell  material,  i.e.,  =  1.5  E^,  we  repeated  the  calculations  and  found  the 

following  results: 

Pj  =  +  0.9842 p  Ib/in. 

/^  =  +  0.8696  p  Ib/in. 

Pj  =  +  0.4166  p  Ib/in. 

Pj  =  -  0.2449  p  Ib/in. 

Pj  =  +  0.02242 p  in. -Ib/in. 

Pj  =  “  0.01315  p  in.-lb/in. 

E^w^  =  +  26.406  p  Ib/in. 

E^w^  =  +  39.609  p  Ib/in. 

=  +  25.919  p  Ib/in. 

E^w^  =  +  38.878  p  Ib/in. 

E^wP  =  +  82.291  p  Ib/in. 

=  -  3.2211p  Ib/in. 

=-  3, 894  p  Ib/in.  2 

max 

ffj  =  -  13.303p  lb/in.2 

*max 


1 

P®* 

ps> 

Outer  Shell 

1  - 10.256  p 

-  9.442  p 

-13.603  p 

-  9.600 p 

Inner  Shell 

1  -  9.565p 

-11.658p 

-12.144p 

- 12.390  p 
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APPENDIX  A 


DERIVATION  OF  THE  FUNCTIONS  At2l,  A^,  A^sl,  AND  A^®! 


If  the  beam-column  effect^  due  to  the  axial  portion  of  the  hydrostatic  pressure  is 
neglected,  then  the  differential  equation  governing  the  axisymmetric  elastic  deformations, 
based  on  small-deflection  theory,  of  a  thin-walled  circular  cylinder  is  given  by:^ 


D 


d^v) 

dx* 


+  E  -  w  =  P.  - 

R2 


[A.l] 


The  homogeneous  form  of  Equation  [A.l]  will  be  used  to  derive  edge  coefficients  for  cylin¬ 
drical  shells  of  short  length  in  which  interaction  effects  between  the  two  ends  of  the  shell 
prevail.  Then  we  have 


d*w  h 
D - +  E  -  u)  =  0 

dx* 


[A.  2] 


The  solution^  of  Equation  [A. 2],  which  solution  describes  the  bending  deformations,  can  be 
written  in  the  form: 


ui^{x)  -  Cj^cos  /3*«cosh  jSa;  +  CjSin  j8**cosh  8* 
+  CjCos  ^X'Sinh  Bx  +  C^sin  jS*»sinh  jS* 


[A.  3] 


and  the  first  three  derivatives  of  [A. 3]  are: 


1 


dwf, 

— —  ^  (C 2  +  C^)cos  6a:-  cosh  jSa:  +  (C^  -  Cj)  sin  /3a:-  cosh  /8a: 
+  (C^  +  6'j)cos  ^a:-sinh  /3a:  +  (C2-C3)sin  jSa:-  sinh  /3a: 


1 

-  -  -  =  C.  cos /8a:- cosh /3a:  -  C,  sin /3a:- cosh /3a: 

2/8  2  dx^ 

+  Cj  cos  /3a:-  sinh  /3a:  -  (7  j  sin  6a:-  sinh  /3a: 


[A.4] 


1  d^Wf, 

-  -  -  =  (Cj  -Cj)  cos  /3®-  cosh  6a:  -  (C^  +  Cj)  sin  /3a:-  cosh  /8a: 

2,8^  dx^ 

+  (C4  -Cj)cos  /3®-sinh  Bx  -  (Cj  +  ^’3)  si"  /3a:-  sinh  j8* 


where  in  Equations  [A.3]  and  [A.4]  we  have 


y3(i-v2) 

/Rh 
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The  integration  constants  ^2’  ^3*  and  appearing  in  Equations  [A. 3]  and  [A. 4] 
will  be  determined  from  a  consideration  of  the  load  boundary  conditions  at  the  edges  of  the 
shell  element;  see  Figure  4.  The  longitudinal  bending  moment  and  the  transverse  shearing 
force  are  related  to  the  derivatives  of  by  the  following  equations: 


dx^ 

dM^  d^w^ 

<2^  - - - 

dx  dx^ 


[A.  5] 


With  reference  to  Figures  4  and  7,  let  it  be  prescribed  that  the  load  boundary  con¬ 
ditions  are  given  by: 


Figure  7  —  Edge  Shears,  Moments, 
Deflections,  and  Rotations  For 
a  Cylindrical-Shell  Element 


at  *  =  +  -^  :  = 

at  a:  =  -  -Q^  = 


[A.6] 

[A.7] 


No  considerations  of  symmetry  with  respect  to  the  point  2  =  0  have  been  taken  in  writing  the 
solution  Equation  [A. 3],  and  in  formulating  the  boundary  conditions.  Equations  [A.6]  and  [A.7]. 
The  development  to  follow  will  be  general  in  this  sense. 

The  substitution  of  Equations  [A.4]  and  [A.5]  into  the  boundary  conditions,  Equations 
[A.6]  and  [A.7],  results  in  the  following  four  equations: 


/3f  .  /3/  pi  pi 

-  =  Cj.  cos  —  •  cosh - Cl  sin  —  •  cosh  — 

21)q2  '*2  2^2  2 

pi  ..  pi  .  pi  ..  pi 

+  C-  cos  —  •  sinh - (7,  sin  —  •  sinh  — 

2  2  2^2  2 
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2D/3’ 


/3Z  Bl  Bl 

=  (Cj  - C3)  cos  —  •  cosh  —  -  (C^  +  Cj) sin  —  •  cosh  — 
2  2  2  2 

Bl  Bl  Bl  Bl 

+  (C4  -  Cj)  cos  —  •  sinh - (Cj  +  C3)  sin  —  •  sinh  — 

2  2  2  2 


Bl  Bl  Bl  Bl 

-  =  C.  cos  —  .  cosh  —  +  C-  sin  —  ■  cosh  — 

20/32  ^2  2^2  2 

^  Bl  ,  Bl  Bl  Bl 

-  C«  cos  —  •  sinh  —  -  C’,  sin  —  •  sinh  — 


[A.8] 


~'^L  Bl  Bl  Bl  Bl 

-  =  (Cj  -  C3)  cos  —  •  cosh  —  +  (C^  +  Cj)  sin  —  •  cosh  — 

2DB^  2  2  2  2 

_  ^  Bl  .  Bl  Bl  Bl 

-  {C ^  ~  cos  —  •  sinh  —  -  (Oj  +  C3)  sin  —  •  sinh  — 

2  2  2  2 


Solving  Equations  [A.8]  simultaneously  gives  the  following  expressions  for  the  four 
integration  constants  Oj,  Cj,  C3,  and  C^: 

(Qr  +  Ql'I  Bl  Bl 

Cj(sinh/3f  +  sin/8f)  =-  -  cos —  ^cosh  — 

2DB^  2  2 

Bl  Bl  Bl\ 

+  -  (cos  —  -  sinh  —  -sin  —  .cosh  —  j 

2p^2  \  2  2  2  2  / 


(Qr-^Ql)  fii  Bl 

(sinh /3f  +  sin /3f)  =-  -  sin  —  •sinh  — 

2DB^  2  2 


(Mr+Ml)  /  Bl 


+  -  (cos  —  -  sinh  -t-  -i-sin  ^  -cosh 

20/3^  \  2 


Bl  Bl 

—  -I-sin  — 
2  2 


E 

2 


[A.9] 


(Qr-Ql)  Bl  Bl 

C«  (sinh  fli-sin  fli)  -  sin  —  ‘Cosh  — 

'  2DB^  2  2 

(Mi-Mr)  .  pi  pi  pi  pi  . 

-  - I  cos  —  •  cosh  —  -I-  sin  —  •  sinh  —  1 

2Dp2  \  2  2  2  2  j 


Bl  Bl 

C3  (sinh /3f- sin /3f)  =-  -  cos  —  -sinh  — 

2pp3  2  2 

(^l-^r)  f  Bl  ^  Bl  Bl  Bi\ 

2pp2  \  2  2  2  2  } 
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If  Equations  [A. 9]  for  the  integration  constants  are  substituted  into  the  deflection 
function,  Equation  [A. 3],  and  the  resulting  expression  is  then  evaluated  at  the  two  edges  of 

the  shell  element,  i.e.,  at  x  =  -  —  for  the  right  and  left  edges,  respectively  (see  Figure  7), 

2 

the  following  equation  is  obtained: 

2  2  r  1  /  (Qr+^l) 

(sinh'^  jSZ-sin"^ ^  + —  =-  -  (cosh  +  cos  iSZ) (sinh  ^f-sin  61) 

2 

(Qr-Ql) 

+  -  (cosh  j8Z- cos  j8Z)(sinhi3Z  + sin  6Z) 

(Mr^M^) 


[A.  10] 


40/32 

(Ml-^r) 

4062 


(sinh  6Z-sin  /3Z)'‘ 


(sinh  6Z  +  sin  /3Z)^ 


The  bending  component  of  the  total  deflection  at  the  right  (R)  edge,  i,e.,  at  a;  =  +  —  ,  is 

2 

then  expressed  in  terms  of  the  applied  edge  shears  and  edge  moments  by  the  following: 


Qf^  /  cosh /3Z  •  sinh /3Z  -  cos /SZ  •  sin /3Z 
20^^  '  sinh2/3Z  -  sin2  6Z 

/  cosh /3Z* sin /3Z  -  sinh /3Z 'COS /3Z 


206^  '  sinh2j8Z  -  sin2  |8Z 
iWyj  /  sinh2 /3Z  +  sin2j3Z  \ 


20,62  \  sinh2  /3Z  -  sin2  jSZ 
/  2sinh  /3Z  •  sin  /3Z 


206^  '  sinh^ /3Z  -  sin^ /3Z 


[A.ll] 
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The  bending  component  of  the  total  deflection  at  the  left  (L)  edge,  i.e.,  at  x=-  —  ,  is  then 
expressed  in  terms  of  the  applied  edge  shears  and  edge  moments  by  the  following: 


Kiz,=- 


Ql  /  cosh  •  sinh /3f  -  cos  jSZ  •  sin 


20/3^ 

Mr 

2D8^ 


sinh^  -  sin^  j3l 

cosh  131 .  sin  31  -  sinh  31  •  cos 
sinh^jSZ  -  sin^)3Z 

sinh^  81  +  sin^  /3Z 


[A.12] 


sinh^  6Z  -  sin^/^Z 

2sinh/3Z  •  sin  81 
sinh^  jQZ  -  sin^  jSZ 


Comparing  the  terms  of  Equations  [A. 11]  and  [A.12]  with  the  corresponding  terms  of 
Equation  [1]  and  the  appropriate  edge  coefficients,  Equations  [7a],  shows  that 


A^^\l3l)  = 


sinh^  81  -  sin^  81 

cosh  /3Z  •  sin  j8Z  -  sinh  pi  •  cos  pi 
sinh^  /JZ  -  sin^  pi 


a[i](^Z)  =  sinh^^Z  +  sin^flZ 
sinh^isz  -  sin^/SZ 


[A.  13] 


a[4]^^Z)=  gsinh  flZ  .  sin/3Z 
sinh^  pi  -  sin^  pi 

The  other  two  lambda  functions,  namely,  A^^]  and  A^^^,  enter  into  the  equations  for  the 
edge  rotations  of  the  shell  element,  and  expressions  for  these  two  functions  are  derived  next. 
When  Equations  [A. 9]  are  substituted  into  the  first  derivative  or  slope  function.  Equation 

[A. 4],  the  slopes  at  the  two  edges  of  the  shell  element,  i.e.,  at  a;  =  ±  —  for  the  right  and  left 
edges,  respectively  (see  Figure  7),  are  given  by:  ^ 
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(sinh^  jSZ-  sin^  ^1) 


2 


(Qr  -Ql) 

40/32 


(sinh/3^  +  sin/3Z)2  + 


{Qr^Ql) 

40/S2 


(sinh  /3Z-sin  + 


20/3 


(cosh  jSZ  +  cos  j8Z)(sinh  ^l  +  sin^l)  + 


[A.14] 


20/3 


(cosh  jSZ-cos  /8Z)(sinh  j3l ~ain  61) 


For  those  terms  in  Equation  [A.14]  which  have  the  double  signs,  it  is  intended  that  the  upper 

sign  apply  to  the  right  edge  and  the  lower  one  to  the  left  edge^. 

Thus  the  rotation  at  the  right  (R)  edge,  i.e.,  at  a;  =  +  —  ,  is  expressed  in  terms  of  the 

2 

applied  edge  shears  and  edge  moments  by  the  following; 


^  sinh2 /31  +  sin2  /31  ^  ^  2  sinh /SI  •  sin  Ql 

20/32  \sinh2/31-sin2/31  f  20/32  V  sinh2  |81-sin2  Ql 

^R  j 

f  cosh  /31  •  sinh  /SI  +  cos  /31  •  sin  /31  \ 

15^  ' 

^  sinh2/3l  -  sin2/31  / 

f  cosh  /31  •  sin  /31  +  sinh  /31  •  cos  /31  V 

D(3  ’ 

^  sinh2  pi  -  sin2  pi  • 

+ 


[A.15] 


and  the  rotation  at  the  left  (L)  edge,  i.e.,  at  ar  - - ,  is  expressed  in  terms  of  the  applied 

2 

edge  shears  and  edge  moments  by  the  following: 


~dw^~ 

1 

!  sinh2  /31  +  sin2  Ql  ^ 

1.  1 

/  2  sinh  pi  •  sin  pi  \ 

dx 

/  =  +  y^-  + 

*  =  --  20^2' 

^  sinh2  /31-sin2  pi  1 

'  20,Q2 

\ sinh2 |81-sin2 j 

M /  cosh /31  •  sinh  +  cos /31- sin  j81 


-L  I 

~  Dfi  \ 


[A.16] 


sinh2/31  -  sin2  /31 
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Mu  /  cosh /3Z»  sin /3Z  +  cos  •  sinh /3Z  \ 

'  sinh^  -  sin^ /Si  ' 

Comparing  the  terms  of  Equations  [A.  15]  and  [A,  16]  with  the  corresponding  terms  of 
Equation  [2]  and  the  appropriate  edge  coefficients,  Equations  [7a],  shows  that  the  remaining 
two  lambda  functions,  besides  those  defined  by  [A. 13],  are  given  by: 


r_i  cosh /Si*  sinh  6i  +  cos  i8i  •  sin/3i 

A*-2J(/Si)  =  - - - ^ - : - — 

sinh^  /Si  -  sin^  Bl 


A^®\/3i) 

sinh^  /Si  -  sin^  /Si 


[A. 17] 


Consequently,  the  set  of  functions  defined  by  Equations  [A.  13]  and  [A.  17]  are  exactly  those 
given  as  Equations  [8]  earlier  in  the  report. 
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APPENDIX  B 


DEVELOPMENT  OF  EDGE  COEFFICIENTS  AND  EXPRESSIONS  FOR  THE 
RADIAL  AND  TANGENTIAL  STRESSES  OF  A  CIRCULAR  ANNULUS 


With  reference  to  Equation  [5]  and  Figure  5,  it  has  been  assumed  that  the  edge  moments 
,  axial  thrusts  and  surface  pressure  p  do  not  give  rise  to  any  radial  displacements 
in  the  plane  of  the  circular  annulus.  Only  the  inplane  forces  and  acting  on  the  outer 
and  inner  circular  boundaries,  respectively,  of  the  annulus  give  rise  to  such  deformations. 

On  Page  418  of  Reference  6  the  following  expression  is  given,  based  on  the  Lame  or 
plane-strain  solution  for  a  thick-walled  tube  subjected  to  simultaneous  internal  pressure 
and  external  pressure  p^: 


M(r)  = 


r2  (1  -  v)  (p .  f  2  +  (1  +  v)  (p,.  -  pp)  r.^r  2 

fP(r2  -  r2) 


[B.l] 


Where  r-  and  are  the  radii  to  the  inside  and  outside  circular  boundaries,  respectively,  of 
the  tube,  and  p^.  and  Pq  are  the  radial  pressures  acting  on  the  inside  and  outside  surfaces, 
respectively,  of  the  tube.  The  variable  “r”  is  the  radial  distance  from  the  axis  of  the  tube 
to  a  point  in  question  through  the  thickness  of  the  tube  wall. 

Adapting  the  solution.  Equation  [B.l],  to  the  present  problem  of  the  circular  annulus, 
we  see  that 


=  P. 


Po  =  + 


[B.2] 


Substituting  [B.2]  into  [B.l]  and  adapting  the  sign  convention  of  Figure  5  for  positive  radial 
displacement,  we  obtain  the  following  results: 


_  f  (1-  y)  {Hf  « .2  +  P  .2)  +  (1  +  v) 

-w(r)= - — - — - 

rE^t(R?  -  P.2) 


[B.3] 


To  find  the  edge  coefficients  gf ,  g.'j ,  g^ ,  and  appearing  in  Equation  [5]  and  its  counter¬ 
part  in  which  i-*j  and  j-*i,  it  is  only  necessary  to  substitute  the  following  successive  four  con¬ 
ditions  into  the  basic  solution,  Equation  [B.3]: 
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To  get  : 

set  T  =  R-; 

=  0 ;  =  1 

[B. 

4] 

■■ 

II 

-*1 

Hf  =  I ;  =  0 

[B. 

5] 

gf  : 

T=Rj\ 

=  1  ;  Hf  =  0 

[B. 

6] 

r^Rj-, 

Hf  =  0  ;  Hf  =  1 

[B. 

7] 

Thus  conditions  [B.4],  [B,5],  [B.6],  and  [B,7]  when  substituted  into  Equation  [B.3]  lead  to 
the  following  equations,  respectively: 


(l  +  u)R^Rj^ 

gA  ^ 

E^t{R?-R?) 

''  E^t(R?-R?) 

(l  +  v)ff.252 

g^  =  - — 

'  E'^t(R^-Rj^) 

il^v)R^R? 

gA  =  - L_ 

E^t{Rf-Rf) 


[B.  8] 


[B.  9] 


[B.IO] 


[B.ll] 


With  i-»l,  ;-»2,  and  t  =  Equations  [B.8],  [B.9],  [B.IO],  and  [B.ll]  become  exactly  Equations 

[10a],  respectively. 

For  expressions  [6]  for  the  edge  rotations  of  the  annulus,  symmetry  and  loading  con¬ 
ditions  for  a  typical  bay  of  the  web-stiffened  sandwich  cylinder  far  removed  from  end  effects 
dictate  that  these  rotations  not  only  total  zero  but  each  and  every  component  is  zero.  The 
more  general  case  shown  in  Figure  5  and  reflected  by  Equations  [6]  will  be  considered  in  a 
separate  report. 

Equations  [18]  and  [19]  for  the  radial  and  tangential  stresses,  respectively,  in  the  cir¬ 
cular  annulus  are  derived  next  by  following  the  solution  given  on  Pages  415  to  418  of 
Reference  6.  The  plane-strain  theory  applied  to  the  axisymmetric  elastic  deflections  of  a 
thick-walled  tube  results  in  the  following  differential  equation: 


1  du 
dr^  T  dr 


[B.12] 
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where  u  is  the  radial  displacement  at  a  point  in  the  wall  of  the  tube  at  a  radial  distance  “r" 
away  from  the  axis  of  the  tube.  The  solution  of  Equation  [B.12]  is  given  by: 


B 

u{t)  =  Ar  +  — 

r 


[B.13] 


The  integration  constants  "A”  and  “fi”  are  determined  from  the  following  deflection  boundary 
conditions: 


at  r  =  R-:  u  =  - 


r  =  R.:  u  =  - 


Substituting  the  conditions  [B.14]  into  the  solution  [B.13]  gives: 


[B.14] 


fi^^Ri-wfRj\ 

\  R?  -  R?  / 


Rf  -  R. 


=  -  H ^ — - 

'  \  R?-R?  ) 


B 


[B.15] 


In  Reference  6  the  radial  and  tangential  stresses,  respectively,  as  a  function  of  the 
distance  “r,”  are  as  follows: 


E^  1 

du 

“1 

II 

1 

__dr 

E^ 

~u 

du"j 

- 

(l-v2) 

rfr  j 

[B.16] 


[B.17] 


Substituting  the  deflection  u  (r)  and  its  first  derivative 
Equations  [B.16]  and  [B.17]  yields: 


a,  =  E^ 


r— 


B 


)  (!  +  •') 

B 


V)  (1  +  v) 


dr 

■i] 

-a 


from  Equation  [B.13]  into 


[B.18] 


[B.19] 
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Note 

4  ^ 

o  +  a.  =  2E  -  =  constant 

which  is  a  consequence  of  the  plane-strain  assumption,  i.e., 


[B.20] 


V 

- - (Pr'^  =  constant 


[B.21] 


Equation  [B.21]  results  by  putting  the  axial  stress  o ^  equal  to  zero  in  the  three-dimensional 
Hooke’s  law. 

The  maximum  radial  stress  occurs  on  the  outer  boundary  of  the  annulus,  i.e.,  at 

T  -  R.^  whereas  the  maximum  tangential  stress  occurs  on  the  inner  boundary,  i.e.,  at 

r  =  Rj.  This  together  with  i-*\  and  j-*2  results  in  Equations  [18]  and  [19]. 
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APPENDIX  C 


DEVELOPMENT  OF  EXPRESSIONS  FOR  THE  SHELL  STRESSES, 
EQUATIONS  [11]  THROUGH  [14] 


In  Reference  7  Salerno  and  Pulos  developed  a  theory  for  the  axisymmetric  elastic 
deformations  and  stresses  in  a  ring-stiffened,  perfectly  circular  cylindrical  shell  subjected  to 
uniform  external  hydrostatic  pressure.  Equations  developed  by  these  authors  for  the  critical 
shell  stresses  are'reviewed  here  and  adapted  to  the  present  problem  of  the  two  coaxial  cyl¬ 
inder  elements  comprising  the  web-stiffened  sandwich  cylinder  structure;  see  Figure  1. 

From  symmetry  considerations  (Figure  4),  the  general  solution  for  the  bending  defor¬ 
mations,  i.e..  Equation  [A. 3],  simplifies  to: 

w^^{x)  =  C ^  cos  ^a;-  cosh  fix  +  C ^  sin  /3a;*  sinh/Sa;  [C.l] 

The  particular  integrals  to  the  differential  Equation  [A.l],  which  constitute  the  membrane 
deformations  and  which  must  be  added  to  the  bending  component  [C.l]  to  get  the  total  deflec¬ 
tion,  are  given  by  Equations  [16]  for  the  outside  and  inside  cylindrical  shells,  respectively. 
The  loading  condition  to  which  Equations  [16]  apply  is  shown  in  Figure  1.  The  total  deflec¬ 
tion  can  thus  be  written  in  the  following  form  to  apply  to  both  cylinders: 

w{x)  =  Wfjix)  +  u)P  =  Cj  cos  .6a;*  cosh  /3a;  +  sin  /3a;*  sinh  jSa;  +  U)P  [C.2] 

The  first  derivative  or  slope  expression  is  then  given  by: 


dw  (x) 
dx 


=  (Cj  +  C^)  /3  cos  /3x*  sinh  jSx  -  (C  j  -  C^)  p  sin  fix  •  cosh  /3x 


[C.3] 


The  integration  constants  Cj  and  are  determined  from  the  following  deformation 
boundary  conditions: 


1 

at  X  =  -  :  w.  =  to/  ;  -  =  +  0/  =  0  [C.4] 

2  dx 

When  the  conditions  [C.4]  are  substituted  into  Equations  [C.2]  and  [C.3],  the  constants  Cj 
and  are  found  to  be: 
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where  the  lambda  functions  and  A^®^  are  defined  by  Equations  [8]. 

The  principal  stresses  in  the  longitudinal  and  circumferential  directions  of  the  shell 
elements  are  given  by  the  following  expressions,  respectively: 


=- 

A 

[C.6] 

w 

P 

[C.7] 

h  ** 

where  the  first  term  in  Equation  [C.6]  and  the  first  two  terms  in  Equation  [C,7]  are  the 
corresponding  membrane  stress  components  and  the  remaining  terms  are  the  bending  components. 
In  terms  of  the  shell  curvatures,  the  stress  expressions  [C.6]  and  [C.7]  become  (see,  for 
example.  Reference  7); 


P  ,  Eh  d^w(x) 

(®)  =  -  -T  -  •  - — 

2(1-v2)  dx^ 


[C.8] 


(®)  —  ^ 


w(x)  P 

-  —  1/  — 

R  h 


^  vEh 

2(1-v.2) 


d^w{x) 

dx^ 


[C.9] 


Substituting  the  deflection  w{x)  and  the  second  derivative  of  w{x)  from  Equation  [C.2]  into 
the  above,  gives  the  following  equations: 
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[-  Cj  sin  (3x  •  sinh  j3x  +  C ^  cos  jSx  •  cosh  /3x] 


[C.IO] 


a^(x)  =  - 


P 

h 


^  Eh^'^ 
(1-^2) 


-E 


cos  ^x  •  cosh  /3a;  + 


sin  jSa;  •  sinh  fix 


[C.ll] 


Once  the  constants  and  as  given  by  Equations  [C.5]  are  substituted  into 
Equations  [C.IO]  and  [C.ll],  the  distributions  of  total  longitudinal  and  total  circumferential 
stress  between  adjacent  supporting  elements,  imposing  the  restraint  conditions  defined  by 
[C.4]  on  the  shell  “edges,”  can  then  be  determined. 

Of  particular  interest  are  the  critical  stresses  that  occur  at  a  point  between  adjacent 
supporting  elements,  i.e. ,  at  a;  =  0,  and  immediately  at  a  supporting  element,  i.e.,  at 

X  =  -  —  ;  see  Figure  4.  These  critical  stresses  are  found  from  Equations  [C.IO]  and  [C.ll] 
to  be; 

AT  MIDBAY  (**0): 


P  +  Eh^‘ 


^  (1-^)  ' 


[C.12] 


=  ~ 


E 

ft 


—  w’’  ~u  —  +  E 


P 

7 


■1  +  yhfi^ 


R 


(l-v2)  J 


[C.13] 


AT  A  SUPPORT  {*  =  ±  -  ): 


P  .  Eh^'‘ 


Ov(= - - - 


^  (l-v2) 


-Cj  si 


sin  —  •  sinh  —  +  C .  cos  —  •  cosh 


] 


[C.14] 
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[C.15] 


E 


-E 


cos  —  • cosh  — 
2  2 


sin  —  •  sinh  — 
2  2 


+ 


The  substitution  of  Equations  [C.5]  for  the  constants  Cj  and  into  Equations  [C.12] 
through  [C.15],  and  the  introduction  of  the  lambda  functions  defined  by  Equations  [8]  into  the 
resulting  expressions  lead  to  Equations  [11]  through  [14],  respectively,  for  the  critical  shell 
stresses. 

I 

The  total  deflections  at  midbay,  i.e.,  a;  =  0,  and  at  a  web  stiffener,  i.e.,  ® 
be  found  by  using  Equations  [C.2]  and  (C.5].  These  are,  respectively: 


to(0)  =  -  (wP  -«)*) 


[C.16] 


(wP-w*) 
\  2  / 


iSZ  /3Z 

cos  —  • cosh  — 
2  2 


+  wP 


[C.17] 
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APPENDIX  D 


DETERMINATION  OF  AXIAL-PRESSURE  LOAD  DISTRIBUTION 
TO  THE  TWO  COAXIAL  CYLINDERS 


With  reference  to  Figure  8,  if  it  is  assumed  that  the  web  elements  do  not  resist  any 
axial  load,  then  force  equilibrium  in  the  longitudinal  direction  requires  that 


Figure  8  -  Distribution  of  Axial-Pressure 
Load  to  the  Two  Cylindrical  Shells 


Pj  Pj  +^2^2 


[D.l] 


where  the  axial  stress  forces  P^  and  in  the  outer  and  inner  cylindrical  shells,  respectively, 
are  the  unknown  quantities  to  be  determined.  To  find  explicit  expressions  for  Pj  and 
another  relationship  between  these  quantities,  the  applied  pressure  p,  and  the  geometry  of 
the  shells  is  needed. 

If  it  is  assumed  that  both  cylindrical  shells  displace  the  same  amount  longitudinally, 

i.e., 

Uj  =  Uj  [D'2] 

then  the  integral  of  the  longitudinal  midthickness  strains  over  a  stiffener  spacing  for  each  of 
the  two  cylindrical  shells  must  be  equal.  Therefore, 


1/2  1/2 

Equation  [D.3]  is  a  consequence  of  the  strain-displacement  relation 


[D.3] 


€ 


X 


du 

dx 


[D.4] 
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Introducing  the  two-dimensional  Hooke’s  law 


1 

[D.5] 

1 

‘0  =  7 


into  Equation  [D.3],  we  obtain 

[(l-v2)a^^-v£:*t^^]^0ar  [D.6] 

Since  the  longitudinal  membrane  stress  and  the  circumferential  membrane  strain  in  each  of 
the  two  shells  are  given,  respectively,  by: 


P 

o.M  =  7 

w 

*  4>^  ~  ~R 


then  Equation  [D.6]  becomes: 


- - -  0*  =  I 

Jo  L  «i  J  Jo  L 


(l-v2)P, 


-  V 


V)^{X) 

li. 


dx 


[D.7] 

[D.8] 


[D.9] 


Substituting  the  deflection  function  [C.2]  together  with  the  appropriate  expressions 
for  the  membrane  deflections  for  each  of  the  two  shells  from  Equations  [16]  into  Equation 
[D.9],  carrying  out  the  indicated  integrations,  and  finally  introducing  Equations  [24],  [22], 
and  [21]  into  the  resulting  expression,  we  obtain 
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Equations  [D.l]  and  [D.IO]  constitute  two  equations  in  the  two  unknown  forces  Pj  and 
/j ;  when  they  are  solved  simultaneously,  the  following  expressions  are  found: 


1  «,Aj  -I  i/£:*A,  vE‘h^ 


[D.ll] 


p  «,  1  -1 

pR^  -v(i-a,)—  - - 

«j  2  Rj  J  ftj 


R,A, 


vE’R^h^ 

Hi 


[D.12] 
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where 


^1  = 


[D.13] 
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